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Abstract 



We reconsider the derivation of the Michael lattice sum rules, which relate the energy and 
action stored in a flux tube of a quark-antiquark pair to the static interquark potential, 
and show that they require essential corrections. We then flnd, using the coupling constant 
sum rule of Karsch, that the total Minkowski fleld energy does not match the interquark 
potential, if one follows conventional notions. The implications of this result are discussed. 



Monte Carlo simulations of non-abelian gauge theories are usually carried out on a 
lattice with equal lattice spacings in the spatial and euclidean time directions. A lattice 
regularization using different cutoffs for the space and time directions must of course 
yield the same results for physical observables in the continuum limit. A knowledge of 
the corresponding regularized action is required to relate thermodynamical observables 
to expectation values of lattice operators, since such a regularization allows one to vary 
independently the temperature and volume of the system [1-3] . The action then depends 
on two coupling constants associated with the temporal and space-like plaquettes. These 
couplings are functions of the spatial lattice cutoff and the anisotropy parameter ^, defined 
as the ratio of the spatial to temporal lattice spacing [2,3]. When taking the continuum 
limit with the anisotropy parameter held fixed, physical observables, such as the interquark 
potential and particle masses, should not depend on ^. By requiring that the string tension 
obtained either from the expectation value of a space-time or a space-like Wilson loop be 
invariant under changes in the anisotropy parameter, Karsch [3] was able to show that 
the sum of the derivatives of the inverse coupling constants squared with respect to the 
anisotropy parameter, evaluated on an isotropic lattice, is given in the continuum limit 
by the first coefl&cient in the perturbative expansion of the /^-function. Furthermore, by 
requiring that in the continuum limit the effective action be independent of the lattice 
regularization chosen, he was able to determine the dependence of the couplings on the 
anisotropy parameter. 

The formulation of SU (N) gauge theories on an anisotropic lattice has been used in 
[4] to derive sum rules relating the potential of a quark- antiquark pair to correlators of the 
action, or of the chromoelectric and chromomagnetic field energy, with the Wilson loop. 
These lattice sum rules are known as Michael's sum rules. Comparing the action sum rule 
given in [4] with the corresponding sum rule in ref. [3], one finds that they disagree by a 
factor of 2. That the action sum rule in [4] is actually incorrect has been noted recently 
in [5], where the authors derive this sum rule within the framework of the continuum 
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formulation. Motivated by these observations, we have reexamined the derivation of the 
Michael sum rules, and find that there are important corrections to both the action and 
the energy sum rule of ref. [4]. Using the coupling constant sum rule of Karsch [3], we 
then find, following conventional lore, that the sum rule relating the interquark potential to 
the field energy in the continuum formulation would be violated in the lattice regularized 
theory. The implications of this result are discussed. 

Consider the ground state energy Eq of a quark-antiquark pair separated by a distance 
R. Quantities denoted with a "hat" will always be understood to be measured in units 
of the lattice spacing. The energy Eq can be calculated from the expectation value of the 
Wilson loop with spatial and temporal extension R and T, respectively: 

Eo = - lim i In < W{R, f) > . (1) 

On an isotropic lattice < W{R, T) > is calculated with the action 

s = p{rr + Vs) (2) 

where, for SU{N), is defined in terms of the bare coupling constant by /3 = and 
VtjVs denote the contributions of the time- like and space-like plaquette variables: 



(3) 

2N' 



2N 

n i>j 

Here i,j label the spatial directions, and U^,y{n) is the lattice variable associated with a 
plaquette located in the /xz^-plane at the lattice site n. 

The lattice energy Eq, defined by (1), is a function of R and /3. Since the self-energy 
contributions to Eq, associated with the quark and antiquark, do not depend on R, they 
can be eliminated by considering the difference Eq{R, (3) — Eq^Rq,/]), where Rq, is some 
reference q'^-separation. Then the subtracted qq potential is given by 

[V{R, P)]subtr = - lim i[ln < W{R, f >]subtr , (4a) 

f ^oo T 



where 

< w{R, T) >= J ^^yy{^,\)e _ ^^^^ 

From here on we will always assume that such a subtraction has been carried out, and 
shall drop the subscript " subtr" for simplicity. Following [4] we now take the derivative of 
(4a) with respect to (5. One then obtains 

= Jim ;f^<Vr + Vs >,,--o, (5a) 

op T^oo T 

where the expectation value < O >qq-o is defined generically by 

<WiR,f)0> _ 

< O >gg-o= I . <0> . (56) 

< W{R, T)> 

Hence the /3-derivative of the potential is expressed in terms of correlators of plaquette 
variables with the Wilson loop. These are the correlators which have been measured in 
Monte Carlo simulations to determine the spatial distribution of the energy density in a 
flux tube connecting a quark and antiquark [6] . 

In the limit T oo, the rhs of (5a) can be further simplified. Since for T — oo the 
expectation value of a plaquette variable is invariant under a "time" translation, we have 
that 

<Va>qq-0 ~ f < V'^ >qq-0, (6) 
T^oo 

where is given by an expression of the form (3), with n running over the lattice sites 
on the fixed time slice. With a Wilson loop extending from n4 = — ^ to n4 = ^ , this time 
slice is conveniently chosen to be the 714 = plane. Then 

^m^h ^ ^ < + >^__^^ aJ2^'l<E\x) + B\x) >,,_o, (7) 

where in the last step we have taken the naive continuum limit. Here — and are 
the Minkowski energy densities of the chromoelectric and chromomagnetic fields expressed 
in terms of the euclidean fields. 



We next use the renormalization group to cast the Ihs of (7) in a form involving the 
potential and its derivative with respect to R. In the limit of vanishing lattice spacing "a" 
we have that 

^vf^Jia)) ^V{R), (8) 

where V{R) is the physical interquark potential, and the behaviour of /3(a) as a function 
of the lattice spacing is given, close to the continuum limit, through the renormalization 
group relation 

a = — e ^^''o . (9) 



Here 6o and bi are given by 

117V , 34 / iV \^ 
The invariance of the Ihs of (8) with regard to changes in the lattice spacing leads to 

sp 9mj)^y^^j^^f,smJi, (11) 



dlna dp dR 
where it is understood that this expression is to be evaluated for i? = ^, and with (3 (a) 
determined by (9). Making use of the relation (11), equation (7) takes the following form 
close to the continuum limit 

In the case of a confining potential, V{R,P) = d-{(3)R, this equation reduces to 

which coincides with that obtained in [3] by making use of the relations (3.7), and the 
equation following it in that reference. The second term appearing on the Ihs of (12), 
which gives rise to the factor of two in (13), has been missed in [4]. Hence the Michael 
action sum rule is incorrect, as was also recently observed in [5], where the authors derive 
this sum rule within the continuum formulation of QCD. 



A second sum rule, relating the interquark potential to the Minkowski field energy, can 
be obtained by requiring that a lattice regularization involving different lattice spacings 
in the temporal and spatial directions should lead to the same physical potential as that 
computed from an isotropic lattice. A similar argument has been used in [3] for the string 
tension computed either from a time- like or a space-like Wilson loop, leading to a coupling 
constant sum rule which will play a key role in our discussion. On an anisotropic lattice the 
action involves two couplings, and /3s [1-3] , associated with the time- like and space- like 
plaquette contributions: 

S = prVr + ^PsVs. (14) 

These couplings are conventionally considered to be functions of the spatial lattice spacing 
a and the anisotropy parameter ^ = a/cr, where a,- is the lattice spacing in the temporal 
direction. They are usually parametrized as follows [2,3] 

/^- = ^r^r\ ^r = ^^e, (15) 

where g1{a, 1) = (7^ (a, 1) = gQ{a), and where the explicit dependence on ^ is chosen in such 
a way that the ^ dependence of 5fo-(a,^),cr = s, r, arises only from quantum corrections. 
In the weak coupling limit /^^-(a, ^) and /^^(a, ^) can be related to the bare coupling j3{a) 
on an isotropic lattice [2,3] by 

i (16) 

^4(a, =/3(a) + 2iVc,(0 + O0-^), 

where the ^-dependence of the functions Co-(0, o" = r, s, have been studied in detail in [3]. 

With the action (14), the lattice potential computed from the expectation value of the 
Wilson loop becomes a function of J2, 13 a, and We now require that in the continuum 
limit a — > 0, Or ^ 0, £ = ala^ fixed, the physical potential, F ~ — V, should not depend 
on the choice of ^, if the couplings /3s and are tuned with a appropriately. Hence for 
a — >^ we have that 

= . (17) 



d 



^F(-,/3,(a,0,/3.(a,0 
a \ a 



Noting that 



dV 



=<V'>qg-Q; cr = T,s, 



(18) 



one finds, upon carrying out the differentiation (17), and then returning to the isotropic 
lattice ^ = 1, that 



ViR,(3) = 



< / r >qq-0 + < / s >qq-0 



(19) 



^=1 



Here V{R,P) is the potential in lattice units computed on an isotropic lattice. Again it 
is understood that this relation only holds for R = R/a,Pa- = Pa{0'iO^ ^^e continuum 
limit. Defining 



1 



± 



^=1 



1=1 



(20) 



expression (19) can be written in the form 



(21) 



where the expectation values are computed on an isotropic lattice. Making use of the 
action sum rule (12) one obtains 



V{Rj)+ri+ 



dlna 
dp 



V{R,P) + R 



dV{R,P) 
dk 



rj-<-V'^ + Vi>gg-o 



(22) 



We emphasize that so far our discussion has not involved any perturbative arguments. We 
now want to interprete the rhs of (22) in terms of a continuum physical observable. It 
is at this (and only this) stage where, following standard lore, we make use of the weak 
coupling approximations (16). From the definition (20) and the relations (16) one finds 
that for a — > 

rj_ ^ P{a) + N[cUl) - c;(l)] P{a) , (23) 

/3^oo 

where c' is the derivative of c^. Hence close to the continuum limit 



V{Rj)+r]+ 



dlna 
dp 



oR 



(24) 



where, as explained before, V'^ (cr = r, s) denotes the contributions to the action of plaque- 
ttes with base on a fixed time slice. In the continuum limit f5{—V'^ + V'g) can be identified 
with the Minkowski field energy measured in lattice units. Thus for a ~ 0, 

/3 < -7^; + >,,-_o~ a^a^^ < -^^(f) + B\x) >,,-_o , (25) 

where {-\E'^) and are the electric and magnetic contributions to the (Minkowski) 
field energy densities expressed in terms of the euclidean fields. Now the field energy, after 
subtracting the self energy contributions, is expected to be related to the potential by 

V{R) = ^jd^x< -E\x) + B\x) >,,-_o • (26) 

Hence on the lattice ^-f ,/3(a)^ should also match the total field energy for a — > 0. 
Thus if (24) and (25) holds, then we are led to conclude that the second term appearing 
on the Ihs of (24) should vanish in this limit. This is the case if either V + = 0, 

^ ' oR 

dlna/d(3 — 0, or rj^ = 0. In the first case V{R,P) has the form V ~ a{P)/R, and hence 
does not allow for a confining potential. The second possibility is clearly excluded for 
SU (N) gauge theories. Hence we conclude that in a confining theory like QCD, r]^ would 
have to vanish for (3 ^ oo. But according to (20) and (15), 




Hence the above line of reasoning would lead to the conclusion that the quantity appearing 
in curly brackets should vanish in the limit (3 ^ oo. This disagrees with the result obtained 
in ref. [3], where this quantity was shown to be given by defined in (10). The author 
was led to this result by requiring the invariance of the string tensions computed from 
space-time and space-like Wilson loops under changes in the anisotropy parameter ^. On 
an isotropic lattice the two string tensions extracted in this way should coincide. By 
making further use of the action sum rule, Karsch was led to the above conclusion. 
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Using the value r]+ = Nbo obtained in [3], the approximation for ?7_ given in (23), and 
the one-loop approximation dlna/d(5 = —477^, one finds that, for a confining potential, 
(24) reduces to 

aR^2P< + V', >qq-Q , (27) 

which, with the identification (25), would violate the energy sum rule (26) by a factor of 
two. We emphasize that the origin of the factor two in (27) is the second term appearing 
on the rhs of (21). With ri- ^ — 2N/gQ, eq. (21) takes the form 

- " - 2N r "1 

V{R,/3) ^ -^[< -V',+V: >,,_o <K + V: >,,_o| . (28) 

The contribution proportional to < V!^ + V'g >qq-o is absent in the expression obtained in 
[4]. The energy sum rule would have the naively expected form, V = < —V'^ + V'g >qq-o, 
if one neglects this apparent perturbative contribution. The action sum rule (13) for a 
confining potential, together with the value rj-^ = Nbo taken from [3] , however tells us that 
this contribution is about |V for a confining potential. 

Summarizing, our above discussion has shown, that if the weak coupling relations (16) 
are used to identify the rhs of (22) with the field energy in the continuum limit, then the 
implementation of the energy sum rule (26) on the lattice demands that 77+ vanishes in this 
limit, which contradicts the coupling constant sum rule obtained in [3]. This latter sum 
rule has however been derived without invoking perturbation theory. We are therefore 
rather tempted to conclude that the weak coupling argument leading to (24) with the 
identification (25) is inadequate, and that the lattice expression for the field energy is that 
defined by the right hand side of (21). This stands in sharp contrast to the standard belief, 
that the field energy is related to < —V!^ + Vg >qq-o, which as we have seen only holds 
for a l/i?-potential. For a linearly rising potential, on the other hand, the contribution 
proportional to < P^+P^ >qq-o is essential, and implies that the rhs of (22) approaches in 
the continuum limit only one half of the field energy stored in the fiux tube connecting the 
quark-antiquark pair, irrespective of the SU{N) gauge group. Clearly the lattice energy 
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sum rule deserves further investigations. In particular it would be of interest to determine 
the ;5-dependence of rj- from Monte Carlo calculations of Wilson loops on anisotropic 
lattices. 
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